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ABSTRACT 









/ A 



An investigation was made of the magnitude of the bearing and dis- 
tance errors which result from two passive directional acoustic sensors 
obtaining a crossfix on a target in a medium with a constant horizontal 
velocity gradient. Results are tabulated in Appendix III. For the con- 
ditions considered, the errors were found to be small in relation to the 
target range but became appreciable in absolute value at very large tar- 
get ranges. However, with a knowledge of the horizontal velocity gradient, 
sensor separation, sensor axis orientation, apparent target range and 
apparent target bearing the errors could be predicted. Further, by 
proper sensor orientation, the errors could be minimized. 
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1. Introduction. 



Straight-line propagation of sound is to be expected only if the 
velocity of propagation is constant throughout the medium. Because of 
the heterogeneous nature of the ocean, this condition is usually not met 
and marked deviations from straight line propogation have been observed 
to occur. 

The principle factors that influence the velocity of sound in sea 
water and therefore result in refraction of the sound ray are temperature, 
salinity, and depth. Variations in salinity are of importance near the 

mouths of large rivers where appreciable amounts of fresh water run into 

the sea, in the vicinity of large ocean currents such as the Gulf Stream, 
and in water close to the surface where rain and evaporation will have a 
maximum effect. However, by far the most important factors influencing 

variations in the speed of sound are differences in temperature and depth 

of the points through which the sound ray passes. 

A considerable amount of work has been done investigating the re- 
fraction of sound rays caused by variations of temperature with depth. 

The aim of these studies has been to increase the detection capability 
of sensors by placing them in such a manner that "shadow zones” are either 
eliminated or caused to be as far from the sensors as possible. Further, 
the studies have investigated ways to make optimum use of sound "chan- 
nels” caused by having a layer of water in which the sound velocity is 
lower than that both above and below it. 

However, once a detection has been made by two directional acoustic 
sensors and the resulting fix plotted, very little has been done to de- 
termine by how much this apparent fix is displaced from the actual target 
location due to horizontal temperature gradients. Although these gradients 
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are normally small in relation to the vertical ones, there are certain 
areas near ocean currents, such as the Gulf Stream, where they do become 
appreciable and thus could possibly cause significant bearing errors in 
crossfixes. Further, as the range at which sensors can detect a target 
is increased through technical improvements and making use of "channel- 
ing”, bearing errors could become significant. 

One investigation of the fixing errors which could result from hori- 
zontal temperature gradients considered the case where there was a sharp 
vertical boundary between two bodies of water having different sound 
propagation velocities. ^6 j The next logical condition to consider is a 
constant velocity gradient. Therefore an investigation of the effect of 
a constant gradient is presented in the ensuing sections of this thesis. 

2 . Theoret ical Background . 

Sound energy takes the form of disturbances of the pressure and den- 
sity of some medium. Therefore, the basic relationships between impressed 
forces and the resulting changes in pressure and density are applicable 
to the theory of sound propagation in water. However, the solution of 
the differential equations applicable to pressure transmission phenomena 
requires a knowledge of certain initial and boundary conditions. This 
knowledge is almost always lacking for sound transmission in sea water. 
Fortunately, it is possible to discuss some of the important characteris- 
tics of sound propagation in sea water without considering the frequency 
and wave structure of the sound. For such discussions the wave can be 
thought of as a beam or ray, and the discussion is concerned only with 
the geometry of the situation. 

The ray concept is based on the assumption that sound energy is 
transmitted along curved paths, called rays, which are straight lines in 
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all parts of the medium where the velocity of sound is a constant, and 
which curve according to certain rules in parts where the velocity of 
sound is changing. However, in making use of the ray path approach to 
the solution of sound propagation problems it must be kept in mind that 
certain assumptions are implicit in this technique. If the conditions 
imposed by these assumptions do not apply, the use of the method is not 
valid. These assumptions are: 

(a) The direction of the ray must not change much over a distance 
of one wave length; that is, the radius of curvature of the rays must 
be large in relation to the wavelength. 

(b) The velocity of sound must not change appreciably over the dis- 
tance of one wavelength, 7 J 

The simplest continuous velocity gradient which can be considered 
is one which is a constant. This particular gradient has the interesting 
property of resulting in a sound ray path which is the arc of a circle. 
This fact is easily shown by a consideration of the following: 

Assume that the velocity of sound, V, at a point (x, y), is a linear 

function of the distance of the point from the origin in the x direction; 

i.e. V = V + gx, where V and V are the velocities of sound at some 
P o P ° 

point P and the origin respectively, and g is the horizontal velocity 
gradient, a constant. Figure 1 illustrates the situation under consider- 
ation. 

In the diagram the velocity gradient is in the horizontal (x) direc- 
tion. The gradient at any point P with coordinates (x, y) is, by 
def init ion, 
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By integrating this expression and specifying that the x distances shall 
be measured from a vertical line through the points where the velocity 
is zero it is possible to avoid the introduction of a constant of inte- 
gration and obtain the expression 



(2-1) 



x 



V 



g 



Then, recalling Snell’s Law 



( 2 - 2 ) 



K 



cos 0 

where 0 is the angle measured from 
velocity line in this case) to the 
Snell’s Law constant for this ray, 
bined by eliminating V to obtain 



the boundary of the medium (an iso- 
ray under consideration and K is the 
equations (2-1) and (2-2) can be com- 



(2-3) x = — §— cos 0 

g 

Differentiating this expression yields 



(2-4) dx = - — sin 0 d 0 

g 

But, from elementary analyticel geometry dx = tan 0 dy (note that in this 
situation angles are measured clockwise from the y axis), thus (2-4) be- 
comes 

(2-5) dy = — cos 0 d 0 

g 

Now integrating (2-5) and specifying that vertical distances, y, 
are to be measured from the horizontal line through the point at which 
0=0 yields 

(2-6) y = - — sin 0 

By squaring (2-3) and (2-6), then adding them, the following result 
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Figure 1. General Ray Path Geometry 




Figure 2* Constant Gradient Ray Path Geometry 



is obtained: 



(2-7) x 2 + y 2 = (JL ) 2 

g 

LT 

This is easily recognized as the equation of a circle of radius . 

g 

The circle is centered at the origin specified in the foregoing paragraphs 
as being the intersection of a horizontal line through the point where 
0=0 and a vertical line through the points where V = 0. 

Since the Snell’s Law constant, K, at any point is the ratio of the 
velocity of sound at that point to the cosine of the angle of inclination 
of the ray at that point, it is readily seen that the radius of the circle 
V 

is actually 2- i.e. the ratio of the velocity of sound at that 

g cos 0p 

point to the gradient times the cosine of the angle of inclination of the 
ray at that point. Thus, if the velocity and angle of inclination of the 
ray are known at some point, the coordinates and angle of inclination of 
the ray are completely specified for the total path of the ray. 

By making use of the fact that a constant velocity gradient results 

in rays that are arcs of circles some useful relationships can be ob- 

tained for tracing the path of a sound ray through the water. Consider 
a circle of radius R with the points (x^, y^) and as F ^g ure 

2. By geometry, 

x^ = R - R cos 0^ = R( 1 - cos 0^) 

x 2 = R - R cos ©2 = R(1 - cos 0 2 ) 

Therefore 

(2-8) ax = x 2 - x^ = R( 1 - cos 9 ^) - R(1 - cos 0-^) 

= R(cos 0^ - cos 

Now, for a linear velocity gradient, g, the velocity, V 2 , at a point, 

(x^, y 2 ) , may be found from the following relationship: 
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V 2 = V x + g^x 



where V is the velocity at the point, (x^, y^) • 



Thus 



(2-9) 



a x = 



V 2 - V 1 



Equating (2-8) and (2-9) yields 
( 2 - 10 ) 

Solving this expression for R, 



R(cos 0 x - cos 0 2 ) = 



V - V 
v 2 1 



( 2 - 11 ) 



R = 



v 2 - V 1 



g(cos 0 ^ - COS 62) 



But, from Snell's Law, 



Vi 

(2-12) cos 0-, = — _ cos 0 . 

1 \ T - 9 



Substituting (2-12) into (2-11) yields 



(2-13) 



R = 



Vo - V n 



V 2 - Vj 



g(y“" cos 0 2 " cos e 2 ) S COS 0 2 ( V^~ 



- 1) 



V 2 - V 1 

v i - -' v 2 - 

g cos e 2 ( — ) 



g cos 0 2 



Similarly, 



(2-14) Ay = Y 2 " Yl = R s ^ n ^2 “ R sin ®i = R ( s ^- n ^ 2 _s ^ n 



R [2 sin I (0 2 - 0^ cos \ (0 2 + ® 1 )| 



But, from (2-8) 



Ax 



cos 0^ - cos ©2 
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Thus, substituting this into (2-14) 



cos 0 2 - cos 0^ 



A x 



>[ 



2 sin i(0 2 



V cos I(0 2 + ^)J 






Ax 



] 





(0 



Simplifying this expression by means of a trigonometric identity yields 



Equation (2-15) is the basic equation used in the model described in 
section 3 which is used to trace the sound rays through the water. 

It should be noted that the minimum radius of the circle, as compu- 
ted from equation (2-13), describing the path of a sound ray is obtained 
when cos 0 2 and the gradient, g, are both a maximum. This means that for 
the model described in section 3 the minimum radius will occur when 0=0° 
and the velocity gradient is 0.0005 meter/sec per meter. This radius will 
be ( 0005)(l) — = 3 x 10^ meters. Since the wavelength of a 200 cycle/sec 
signal (approximately the frequency of a low frequency noise source on a 
medium speed submarine [ 4 ]) is about 7.5 meters it is readily seen that 
the conditions of assumption (a), that the radius of curvature must be 
large in relation to the wavelength, are met and the use of the ray path 
method is valid. Further, the conditions of assumption (b) , that the 
velocity of sound does not change appreciably over one wavelength, are 
fulfilled since the velocity only changes by , x .0005 _ 2,5 x io ~6 - 



(2-15) 



y = Ax cot ( — ) 




1500 
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. 00025 % in a distance of one wavelength. 



3. The Model . 

In modeling the problem under investigation the following assumptions 
were made: 

a. A detection had been made by two directional acoustic sensors and 
the apparent fix plotted. 

b. There was no change in horizontal gradient with depth, therefore 
the amount of bending of the ray was the same at any depth and the prob- 
lem could be reduced to one of finding the ray path in a horizontal plane 
only. (It is realized that the horizontal gradient does, in actuality, 
change slightly with depth, thus all results are approximations. How- 
ever, to treat the problem in three dimensions would require an amount of 
time exceeding that available for the conduct of the present investiga- 
tion.) 

c. The velocity gradient was a constant in one direction and zero 
in a direction perpendicular to this one. This allowed the use of the 
circle approximation to the path of the sound rays. 

d. The directional acoustic sensor had no bearing error. To add 
this to the problem would obscure the effect of the temperature gradient. 

Figure 3 is an illustration of the situation studied. An apparent 
target was located at XAPP, YAPP, the apparent fix obtained from sensors 
A and B located at XI, Y1 , and XXI, YY1 respectively. The bearings of 
the apparent target from sensors A and B were BRGA and BRGB respectively. 
To relate the fix to some common reference point, an "origin" was chosen 
at the midpoint of a line joining the two sensors. The bearing and range 
of the apparent target from this point were BRGP and RAPP. However, due 
to the horizontal velocity gradient, some bending of the sound rays 
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Figure 3* Computer Program Geometry 
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Figure 4. Computer Program Geometry (Sensor Axis Rotated) 
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occurred and the actual target was located at XACT, YACT. The bearing 
and range of the actual target from the origin were BRGG and RACT. The 
distance from the apparent target to the actual target was the ERROR and 
the bearing of the actual target from the apparent target was called BRGT. 
RNGERR was the difference between the actual and apparent ranges. 

If the sensors had always been located along an isovelocity line, 
Figure 3, with the bearings and distances shown, would be sufficient to 
illustrate the geometry of the problem. However, the inclination, by an 
angle OMEGA, of an axis through the two sensors with respect to an iso- 
velocity line resulted in a need to define angles with respect to an 
isovelocity line. By selecting the x and y axes perpendicular to and 
along an isovelocity line respectively it was possible to make the best 
use of the fact that the velocity gradient was a constant in one direc- 
tion and zero in the other. These angles measured with respect to the 
isovelocity line were assigned Greek names corresponding to the letter 
appearing after the prefix BRG; i.e. BRGA, the bearing with respect to an 
axis through the two sensors, became ALPHA, the bearing with respect to 
an isovelocity line. The Greek letter bearings equal OMEGA plus the 
appropriate BRG. Figure 4 illustrates the problem when the sensor axis 
is rotated with respect to the isovelocity line. 

Making use of the relationships developed in section 2 relating the 
sound ray path to the arc of a circle, a computer program was written to 
trace the actual path of the sound ray through the water. Thus, given 
initial bearings of the apparent target from two sensors, it calculated 
the actual target position from the point of intersection of the actual 
ray paths. The program is included in Appendix I. 

The "mechanics" of the program are illustrated in Figure 5, a flow 



18 




BRGP = Apparent Target Bearing 

RAPP = Apparent Target Range 

RAPPMAX = Maximum RAPP to be Considered 

SS = Sensor Separation 



Figure 5. Flow Chart for Program 



19 



chart for the program. For the reader interested in greater detail a more 
detailed block diagram may be found in Appendix II. The Roman numerals 
in the blocks of Appendix II correspond to the numerals in the left hand 
margin of the program in Appendix I. 

For a given velocity gradient and sensor separation an apparent tar- 
get position was assumed. From this position the apparent bearings and 
distance of the target from the origin and the sensors was calculated. 

Then the actual position of the target was determined by calculating the 
actual ray paths through the water and computing their point of intersec- 
tion. Finally, the errors were calculated. A new apparent target posi- 
tion was then assumed and the calculations repeated for this new position. 
Apparent target positions were assumed for each 10 degrees of bearing 
from an axis through the sensors (except for the 0 and 180 degree posi- 
tions where no fix is possible) and for various apparent ranges to a 
maximum range of five times the sensor separation. Further, as stated 
earlier, the angle of inclination of the axis through the sensors to the 
isovelocity line was varied in 10 degree increments from 0 to 90 degrees 
for selected sensor separation and apparent range combinations. 

4. Results . 

To evolve a means for estimating the crossfix errors resulting from 
horizontal temperature gradients, the effect of the following parameters 
on the errors was studied using the model described in section 3: 

(a) target bearing 

(b) magnitude of velocity gradient 

(c) sensor separation 

(d) target apparent range 

(e) sensor axis orientation 
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The range of values over which the parameters were varied was as 
follows : 

(a) target bearing (BRGP) : 10° increments from 10° to 350° relative 
to the sensor axis except for 180° where no fix was possible due to target 
alignment with the sensor axis. 

b) velocity gradient (VGRAD) : .0001 to .0005 meter/second per meter. 
These were found to be the approximate maximum velocity gradients which 
exist over long (about 100 km.) distances in the deep sound channel and 
over short (about 40 km.) distances near the surface at the Gulf Stream. 

[*] 

(c) sensor separation (SS) : 1 km. to 400 km. These were considered 
to be the minimum and maximum feasible for obtaining a crossfix. 

(d) target apparent range (RAPP): 1 km. to 2,000 km. but never less 
than .1 nor more than 5 times the sensor separation since beyond these 
values the target bearings from the sensors become almost equal. 

(e) sensor axis orientation (OMEGA): sensor axis 0° to 90° from the 
isovelocity line in 10° increments. 

As might be expected, the errors were found to be small in relation to 
the target distance from the sensors but to be appreciable in absolute 
value for targets at long ranges from the sensors or in areas where a 
large velocity gradient existed. It was found that the errors were pre- 
dictable and that they could be minimized by proper sensor orientation. 

The results are summarized in sections 4.1 through 4.8 and depicted 
graphically in Figures 6 through 20. Although all of the computer out- 
put data used in plotting the graphs for the figures is not included, 
samples of the data which illustrate particular effects are included in 
Appendix III. Also, it should be noted that the abbreviations used in 
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the figures are defined at the beginning of the computer program in Appen- 
dix I. 

4. 1 Effect of Target Bearing . 

Figure 6 illustrates the influence of target bearing (BRGP) on the 
magnitude of the crossfix errors. As the figure and the data in Tables 
1, 2 and 3 show, for a given velocity gradient condition the error was 
very much influenced by target bearing from the sensor axis when the sen- 
sor axis was oriented along an isovelocity line. However, when the sen- 
sor axis was at 90° to the isovelocity line, i.e. along the velocity 
gradient line where OMEGA = 90, target bearing had almost no effect on the 
error. The influence of sensor axis orientation is discussed further in 
section 4.5. For the case where the sensor axis was along an isovelocity 
line the errors were seen to increase rapidly as the target bearings ap- 
proached the sensor axis. As discussed in section 2, this was due to the 
fact that the radius of curvature of the sound ray approaches its minimum 
value as the angle which the ray makes with the isovelocity line approaches 
zero. 

4.2 Effect of Velocity Gradient . 

Figure 6 also illustrates the influence of velocity gradient (VGRAD) 
on the magnitude of the cross-fix errors. As the figure and the data in 
Tables 1, 2 and 3 show, the errors were a linear function of the velocity 
gradient, i.e. the velocity gradient of .0005 m/sec per meter resulted in 
errors 5 times larger than those for a gradient of .0001. Although Fig- 
ure 6 depicts the particular case of a 2 km. sensor separation and 10 km. 
target range, the effect which it illustrates was found to be the same 
for all sensor separation and target range combinations examined. The 
plot for the 90° rotation of the sensor axis from an isovelocity line 
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0 90 180 270 360 

OMEGA. = 90 



VGRAD = .0005 Reference; 

VGRAD = .00025 Tables 1,2,3 

VGRAD = .0001 

SS = 2 km, RAPP = 10 km. 



Figure 6. Effect of Velocity Gradient 






BRGP 
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illustrates the linear relationship between velocity gradient and cross- 
fix errors even more clearly than the plot for the 0° rotation of the 
sensor axis. 

4. 3 Effect of Sensor Separation . 

The effect of sensor separation (SS) on crossfix errors is not as 
simply described as the effect of velocity and target bearing. As shown 
in Figures 7 and 8 and by the data in Tables 4, 5, 6, 7, and 8 it was 
found that for any given target range, sensor separation had much less 
effect on the errors than target bearing. Also, the minimum error due to 
target bearing generally increased and the maximum error generally de- 
creased as the sensor separation was increased. However, when the sensor 
separation was twice the apparent target range the maximum error due to 
target bearing was at a minimum value and the minimum error was at a max- 
imum value, thus there was the least dependence of error on target bear- 
ing. When the sensor separation was increased to values more than twice 
the apparent target range the minimum error continued to increase and the 
maximum error now increased also. The effect of sensor separation on 
crossfix errors when the sensor separation was greater than twice the 
target range is discussed further in section 4.8. 

An examination of the data in the figures and the tables revealed 
no simple relationship between the magnitude of the error and the sensor 
separation per se. However, as will be discussed in section 4.6, it was 
possible to formulate an estimating relationship for errors based on the 
ratio of target range to sensor separation. Also, the large variation in 
error for the case where the sensors were 10 km. apart and oriented 90° 
to the isovelocity line will be discussed in section 4.8. 
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ERROR - METERS 




OMEGA = 0 




. . — ss = 1 km. Reference: 

SS = 5 km. Tables 4,5,6 

SS = 10 km. 

RAPP = 5 km. , VGRAD = .0005 



Figure 7. Effect of Sensor Separation 
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ERROR - KILOMETERS 




SS = 100 km. Reference: 

— - — SS ■ 400 km. Tables 7,8 

RAPP = 500 km., VGRAD = .0001 



Figure 8. Effect of Sensor Separation 

f\ 
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4.4 Effect of Target Range . 



Figures 9 and 10 show the dependence of crossfix errors on target 
range (RAPP) . Samples of the data used to plot the graphs may be found 
in Tables 4, 8-12, and 13. Figure 9a shows the absolute error as a func- 
tion of target bearing for three target ranges and a sensor separation of 
1 km. Figure 9b shows the relative error (the ratio of the error to the 
target range) for the same target and sensor positions. Figure 10 shows 
the relative error as a function of target bearing for different sensor 
separation and velocity gradient conditions. 

In examining the data it may be noted that there was no simple 
linear relationship which related the absolute error at one target range 
to another, thus interpolation is difficult. However, by plotting rela- 
tive error versus target bearing, as in Figure 9b, it may be seen that 
the relative error was a linear function of target range. 

As Figure 10 illustrates, when the sensors were on an isovelocity 
line, this relationship was found to be generally true for extreme sen- 
sor separations and target ranges also. However, for target ranges less 
than one-half the sensor separation and for target bearings less than 45° 
from the sensor axis the relationship did not hold. As was mentioned in 
section 4.3, the case when the target was less than one-half the sensor 
separation will be discussed in section 4.8. At target bearings less than 
45° from the sensor axis there was less than a 10° difference between the 
target bearings from the two sensors, thus it is unlikely that a "fix" 
would be considered from these two bearings due to the large errors which 
would be induced by small differences in the measurement of either bearing. 
When the sensor axis was at 90° to the isovelocity line, the relative er- 
ror was a linear function of target range for all target bearings, but as 
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O 90 180 270 360 BRGP 

Figure 9a. Effect of Target Range on 

Absolute Error (Short Range) 




Figure 9b. Effect of Target Range on 

Relative Error (Short Range) 



RAPP = 5 km. References: 

RAPP = 3 km. Tables 4,9,10 

. RAPP = 1 km. 

SS = 1 km., VGRAD = .0005, OMEGA = 0 
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in the case when the sensor axis was along an isovelocity line, when the 
target range was less than one-half the sensor separation the linear re- 
lationship did not exist. 

4. 5 Effect of Sensor Axis Orientation . 

As was noted in Figures 7, 8 and 10, the orientation (OMEGA) of the 
sensor axis with respect to the velocity gradient line had a very large 
effect on the magnitude of the errors. In general, errors were minimized 
when the sensor axis was aligned with the velocity gradient line. Fur- 
ther, the errors were practically constant for all target bearings when 
the sensor axis was aligned with the velocity gradient line. The only 
exceptions to this occurred when the target range was less than one-half 
the sensor separation and this case is discussed in section 4.8. 

Figures 11, 12 and 13 illustrate the effect of rotation of the sen- 
sor axis from alignment with an isovelocity line (OMEGA equals 0°) to 
alignment with the velocity gradient line (OMEGA equals 9 (P) . The data 
used to construct the graphs may be found in Tables 14 through 23 of Ap- 
pendix III. 

Figure 11 is a histogram showing the maximum and minimum errors, de- 
pending upon target bearing, for targets at an apparent range of 6 km. 
from sensors 2 km. apart. As the graph shows, by rotating the sensor 
axis from 0° to 90° the maximum errors were reduced by a factor of 10 to 
a value approximately equal to the minimum error. Figures 12a and 12b 
further illustrate the changes in the errors brought about by rotating 
the sensor axis. 

Figure 13 depicts the actual and apparent target locations when the 
sensor axis was aligned with an isovelocity line and when the sensor axis 
was aligned with the velocity gradient line. It should be noted that the 
graph is not to an absolute scale, since on such a scale even the maximum 
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Figure 11. Effect of Sensor Axis Orientation 
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Figure 12a. Effect of Rotation of Sensor Axis 
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ERROR - METERS 
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Figure 12b. Effect of Rotation of Sensor Axis 
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error of 68 m. in 6000 m. could not be seen. However, the relative magni- 
tudes of the errors are correct. As Figure 13 and the data in Tables 14 
through 23 show, the actual target was closer than the apparent target 
when the sensors were in a lower velocity medium than the target and fur- 
ther away when the sensors were in a higher velocity medium. When the 
sensors were oriented along the velocity gradient line, however, this re- 
lationship was reversed. 

4. 6 Comparison at Gonstant Reduced Range . 

Reduced range (REDRNG) was defined to be the ratio of the apparent 
target range to the sensor separation. This term was introduced to facil- 
itate estimation of errors for any sensor separation and target range when 
the errors at a different sensor separation and target range are known. 
Figure 14 illustrates the linear relationship between relative errors at 
one sensor separation and those at another when compared at the same re- 
duced range. 

Figure 15 depicts relative error as a function of reduced range for 
various target bearings and a velocity gradient of .0005 m/sec per meter. 
As an example of the usefulness of such a plot, consider the following: 

For a target located at a range of 5 km. on a 30° bearing from the 
midpoint of a line connecting two sensors 1 km. apart, Figure 15 shows a 
relative error of about .003 (REDRNG = 5., BRGP = 30). For a target lo- 
cated at a range of 500 km. on a 30° bearing from the midpoint of a line 
connecting two sensors 100 km. apart (at a reduced range of 5.) a predic- 
ted relative error would therefore be .3 (.003 times 100) for a gradient 
of .0005. Thus, using the fact (from section 4.2) that relative errors 
vary linearly with velocity gradient, for a gradient of .0001 a predicted 
relative error would be .06. The absolute error would then be estimated 
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Figure 13. Actual Target Locations 
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to be 30 km. (.06 times 500). This compares favorably with the actual 
error of 28,151 m. found in Table 7 for this situation. 

Over a wide range of sensor separations, target ranges and velocity 
gradients considered the error estimates obtained in this manner were al- 
ways within 14% of the value computed by the model described in section 3. 
The unusual behavior of the curves in Figure 15 below reduced ranges of .5 
will be discussed in section 4.8. 

4.7 Effect on Bearing Error . 

For targets located such that there was at least a 10° bearing differ 
ence between the bearings from the two sensors, the maximum error in bear- 
ing from the origin on the sensor axis was 2.2° (Table 13). Except for 
the extreme sensor separations and target ranges these bearing errors were 
less than 1, thus a graphic portrayal of these would not be meaningful. 
However, it was noted that the bearing of the actual target from the appar 
ent target position was a function of target bearing from the origin and 
was independent of velocity gradient, sensor separation and apparent tar- 
get range except when the range was less than the sensor separation. This 
case will be considered in section 4.8. 

Figure 16 illustrates the relationship between the bearing of the 
actual target from the apparent target (BRGT) and the bearing of the ap- 
parent target from the origin (BRGP) . Four figures have been drawn to 
illustrate the change in this relationship as the sensor axis was rotated 
from alignment with an isovelocity line to alignment with the velocity 
gradient line. As the figures show, the actual target was always toward 
the lower velocity medium from the apparent target position except for a 
sensor axis orientation along the velocity gradient line. 
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Figure 16a. Target Bearing Corrections 
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BRGT - DEGREES 





Figure 16b. Target Bearing Corrections 
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4.8 Effects When Target Range is Less Than Sensor Separation . 



In the foregoing sections numerous references have been made to the 
fact that when the target range was less than the sensor separation, and 
particularly when the target range was less than one-half the sensor separa- 
tion, many of the estimating relationships were no longer valid. 

As may be seen from Figures 7 and 15, when the sensor axis was aligned 
with an isovelocity line, the effect of target bearing on the distance er- 
rors was the same for reduced range values less than 0.5 as it was for 
those above 0.5. However, when the sensor axis was rotated to alignment 
with the velocity gradient line, instead of the errors being almost con- 
stant the errors were greatest at those target bearings where they had been 
the least when the sensor axis was along an isovelocity line. This may be 
seen in Figure 7 and by a comparison of the data in Table 25 with that in 
Table 26. 

Figures 17a and 17b illustrate the change in the absolute errors as 
the sensor axis was rotated. As may be seen in Figure 17a, the absolute 
errors decreased as the axis was rotated but when the reduced range equaled 
0.5, the variation in error due to target bearing remained constant. This 
is contrasted with the behavior, as indicated by Figures 11 and 17b, of 
the error for values of reduced range both above and below 0.5 where the 
variation in the error as well as its absolute value decreased when the 
sensor axis was rotated to alignment with the velocity gradient line. 

Figure 18 depicts absolute distance errors as a function of reduced 
range for reduced ranges from 0.1 to 1. As the figure shows, the maximum 
errors increased in value as the reduced range increased or decreased from 
a value of 0.5 when the sensor axis was along an isovelocity line but de- 
creased uniformly as the reduced range decreased when the sensor axis was 
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SS = 400 km., RAPP = 200 km. Reference: 
VGRAD = .0001 Tables 25,26 



Figure 17a. Effect of Sensor Axis Orientation 
(Reduced Range = .5) 




SS «= 400 km., RAPP “ 80 km. Reference: 

VGRAD = .0001 Tables 27,28 

Figure 17b. Effect of Sensor Axis Orientation 
(Reduced Range = .2) 
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Figure 18. Effects at Reduced Ranges Less Than 1.0 
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along the velocity gradient line. 

Figures 19 and 20 illustrate the effect of apparent target bearing 
from the origin (BRGP) on the bearing of the actual target from the appar- 
ent target (BRGT) . The figures also illustrate how the bearing correction 
varied with the sensor axis orientation. A comparison of these figures 
and the data in Tables 24-28 with Figure 16 and the data in Table 8 reveals 
that as the reduced range was lowered from a value greater than one (Fig- 
ure 16) through 0.5 (Figure 17) to 0.3 (Figure 18) the relationship between 
the bearing from the apparent to the actual target and the bearing from the 
origin to the apparent target reversed. 

5 . Conclusions and Recommendations . 

As a result of this investigation the following conclusions may be 
drawn regarding the crossfix errors of two directional acoustic sensors 
in a medium with a constant horizontal velocity gradient: 

a. Although the distance errors are small relative to the target 
range, they can become appreciable in absolute value at very large target 
ranges . 

b. For gradients within the values studied, the distance errors are 
a linear function of the velocity gradient. 

c. The distance errors are sensitive to the ratio of the apparent 
target range to the sensor separation. If two directional sensors are 
used to obtain crossfixes and it is known from past experience that a 
large proportion of their detections occur at a certain range, they can 
be placed so as to minimize the crossfix error. 

d. The distance errors are sensitive to target bearing and sensor 
axis orientation. By placing the sensors along the velocity gradient line, 
distance errors due to target bearing can be minimized. 
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e. The distance errors are predictable, i.e. the errors can be esti- 
mated from a knowledge of the horizontal velocity gradient, the sensor 
separation, the sensor axis orientation, the apparent target range and the 
apparent target bearing. 

f. The bearing from the apparent to the actual target is independent 
of velocity gradient, sensor separation, and apparent target range. It is 
only a function of target bearing and sensor axis orientation, thus it is 
also predictable. 

Since the model used to formulate these conclusions did not take into 
account the variations in the path of the sound ray due to vertical velocity 
gradients or horizontal gradients other than linear ones, the investiga- 
tion of errors using a more sophisticated three-dimensional model is recom- 
mended. Such an investigation should, of course, determine the sensitivity 
of the errors to bottom conditions such as slope, since a sloping bottom 
will also result in deflection of the sound ray from its original course. 
Further, since some directional sensor systems depend upon the time phase 
relationship between a pulse received by two sensors, the errors caused by 
the actual paths being of different lengths than the apparent paths may be 
entirely different from those calculated by the model used in this investi- 
gation. Therefore, an investigation of this problem should be considered. 



45 



BRGT - DEGREES 




OMEGA = 0 




0 90 180 270 360 BRGP 

OMEGA = 30 

SS = 400 km. Reference: 

RAPP = 200 km. Table 25 

VGRAD = .0001 •’ 



Figure 19a. Target Bearing Corrections 
(Reduced Range = .5) 
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Figure 19b. Target Bearing Corrections 
(Reduced Range = .5) 
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Figure 20a. Target Bearing Corrections 
(Reduced Range = .2) 
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Figure 20b. Target Bearing Corrections 
(Reduced Range = .2) 
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-COOP * » A DORSEY, S/ 1 S/2S , 3 , 5 000 , 2 . 

-FTN »L » E * 

PROGRAM HRZGRD1 

C THIS IS A PROGRAM TO COMPUTE THE BEARING AND DISTANCE 
C ERRORS IN A SONAR FIX WHICH RESULT FROM HORIZONTAL 
C TEMPERATURE GRADIENTS IN THE WATER. 

C THIS PROGRAM IS WRITTEN IN FTN-63 FOR A CDC 1604 COMPUTER. 

C THE Y AXIS IS AN ISOVELOCITY LINE WHERE THE VELOCITY IS VO. 
C THE ORIGIN (X=0,Y=0) IS A POINT MIDWAY BETWEEN THE TWO 
C SENSORS AND IS THE INTERSECTION OF THE Y AXIS AND A LINE 
C THROUGH THE TWO SENSORS. 

C SENSOR A IS LOCATED AT X1*Y1 AND SENSOR B IS AT XXl,YYl. 

C ALL BEARINGS WITH THE PREFIX BRG ARE RELATIVE TO AN AXIS 
C THROUGH THE TWO SENSORS. 

C ALL BEARINGS DESIGNATED BY GREEK LFTTERS ARE RELATIVE TO 
C THE Y AXIS. 

C ALL RANGES AND DISTANCES ARE IN METERS. 

C ANGLES IS A SUBROUTINE TO CALCULATE ANGLES FROM 0 TO 360 
C DEGREES USING TANGENTS. 

C BRGA IS THE BEARING FROM A TO THE APPARENT TARGET. 

C BRGB IS THE BEARING FROM B TO THE APPARENT TARGET. 

C BRGG IS THE BEARING FROM THE ORIGIN TO THE ACTUAL TARGET. 

C BRGP IS THE BEARING FROM THE ORIGIN TO THE APPARENT TARGET 

C BRGT IS THE BEARING FROM THE APPARENT TO THE ACTUAL TARGET. 
C ERROR IS THE DISTANCE BETWEEN APPARENT AND ACTUAL POSITION. 
C RACT IS THE ACTUAL RANGE OF THE TARGET FROM THE ORIGIN. 

C RAPP IS THE APPARENT RANGE OF THE TARGET FROM THE ORIGIN. 

C RAPPINC IS THE RANGE INCREMENT. 

C RAPPMAX IS THE MAXIMUM APPARENT RANGE CONSIDERED. 

C REDRNG IS apparent RANGE DIVIDED BY SENSOR SEPARATION. 

C RELERR IS THE RELATIVE DISTANCE ERROR ( ERROR /RAPP ) . 

C RNGERR IS THE DIFFERENCE BETWEEN RACT AND RAPP. 

C RPD CONVERTS ANGLES IN DEGREES TO ANGLES IN RADIANS. 

C (RADIANS=DEGREES TIMES .01745329252 RADIANS PER DEGREE) 

C SS IS SENSOR SEPARATION. 

C VO IS THE VELOCITY AT X=0» VGRAD IS THE VELOCITY GRADIENT. 

C XACT »YACT IS ACTUAL TARGET POSITION. 

C XAPP » YAPP IS APPARENT TARGET POSITION. 

C XINC IS THE DISTANCE INCREMENT IN THE X DIRECTION. 

C YDIFF IS THE DESIRED PRECISION IN CLOSENESS OF THE 
C COMPUTED AND ACTUAL POINTS OF. INTERSECTION OF TWO RAYS. 

C ALPHA IS THE BEARING OF THE APPARENT TARGET FROM SENSOR A. 

C BETA IS THE BEARING OF THE APPARENT TARGET FROM SENSOR B. 

C DELTA IS THE DIFFERENCE BETWEEN GAMMA AND PHI. 

C GAMMA IS THE BEARING FROM THE ORIGIN TO THE ACTUAL TARGET. 

C OMEGA IS THE ANGLE BETWEEN THE Y AXIS AND A LINE THROUGH 
C THE TWO SENSORS. IT CAN VARY FROM 0 TO +90 DEGREES. 

C PHI IS THE BEARING FROM THE ORIGIN TO THE APPARENT TARGET. 

C PI IS 3.1415926536 

C TAU IS THE BEARING FROM THE APPARENT TO THE ACTUAL TARGET. 
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C 

C 

c 



I 

II 
III 



IV 




V 



V 



AC 



c 



c 



VI 



c 




READ IN DATA. 

100 RPD = .01745329252 $ P I =3 . 141 5926 536 

110 READ 111, VO ♦ VGRAD » YD I FF 

111 FORMAT ( 3F 10 . 0 ) 

120 READ 121, SS*RAPP*RAPPMAX»RAPPINC»XINC 

121 FORMAT (5F10.0) 

PROGRAM STOPS WHEN SENSOR SEPARATION IS NEGATIVE. 

122 IF (SS) 900,900,125 
125 OMEGA=0 . 0 

132 OMEGAR=OMEGA*RPD 

140 Y1 =SS/ 2*COSF ( OMEGA R ) $ YYl=-Yl 

141 X1=SS/2#SINF( OMEGA R ) $ XXl=-Xl 

170 REDRNG=RAPP/SS 

183 PRINT 184, OMEGA, SS,RAPP*REDRNG, VO, VGRAD, YDIFF 

184 FORMAT (1H1»16(/)»5X» 8 H OMEGA = F3.0,9X»5HSS = F6.0,9X, 
17HRAPP = F7.0,9X,9HREDRNG = F4 . 3 , / / 1 5X , 5HV0 = F6.1, 

2 12X , 8HVGRAD = F7 . 5 , 1 2X , 8HYD I FF = F3.1,//) 

185 PRINT 186 

186 FORMAT ( 2X,4HBRGP4X*4HRACT6X,2HV26X*4HBRGA4X,4HBRGB5X, 
15HERR0R4X,4HBRGT3X,5HDELTA3X,6HRNGERR3X,6HRELERR ) 

200 DO 600 1=1,35 

201 I F ( 1-18) 205,520,205 

ROUTINE TO SET APPARENT TARGET POSITION 

205 BRGP- = 10. *1 

206 I COUNT = 0 / 

210 PHI=OMEGA+BRGP 

TEST TO ENSURE THAT PHI IS LESS THAN 360 DEG. 

211 IF(PHI-360. >213,213,212 

212 PH I =PH I -360 • 

ROUTINE TO COMPUTE RELATIVE BEARING OF TARGET FROM SENSORS 

213 PH I R = PH I *RPD 

214 XAPP=RAPP*SlNF( PHIR) $ YAPP=RAPP*COSF ( PHI R ) 

220 K= 1 $ CALL ANGLES ( X 1 , Y 1 , XAPP , YAPP , ALPHA , K ) 

221 ALPHAR=ALPHA*RPD $BRGA=ALPHA-OMEGA 

222 IF(BRGA) 223,225,225 

223 BRGA=8RGA+360 • 

225 K=2 $ CALL ANGLES ( XX 1 ,YY 1 , XAPP , YAPP , BETA , K ) 

226 BETAR=BETA*RPD $BRGB= BETA-OMEGA 

227 IF(BRGB) 228,230,230 

228 BRGB=BRGB+360. 

ROUTINE TO TRACE RAYS AND COMPUTE ACTUAL TARGET POSITION 

230 XA=X1 $ YA=Y1 . $ XB=XX1 $ YB=YY 1 

231 I F ( OMEGA ) 232,299,232 

232 VA=V0+VGRAD*X1 $ V8= V0+VGRAD*XX1 

ALPHAR1 = ALPHAR $ BETAR1 = BET AR 

TEST TO SEE IF XI IS LESS THAN XAPP IN ABSOLUTE VALUE 

233 I F ( ABSF ( XAPP ) — ABSF ( X 1 ) ) 240*280,280 
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ROUTINE TO COMPUTE XACT ♦ YACT WHEN Xl IS GREATER THAN XAPP 
IN ABSOLUTE VALUE. 

240 RATIOA=COSF( ALPHAR 1 ) /VA S RAT I OB = COSF ( BET AR 1 ) /VB 

X2“XX1 $ XX I NC*X INC 

241 X2=X2+XXINC $ V2=V0+VGR AD*X 2 

I F ( ABSF(V2*RATIOA) .GT. 1. ) 257,242 

242 IF(ABSF(V2*RATIOB) .GT.l. > 257,243 

243 ALPHAR 2 = 2 . *P I-ACOSF ( V 2*RAT I OA ) 

BETAR2=ACOSF (V2*RATIOB ) 

DXA = XA-X2 $ DXB = X 2-XB 

NOTE SIGN OF DX A IS ALWAYS +. 

DYA=DXA/TANF ( ( ALPHAR 1+ALPHAR2 ) /2. ) 

DYB=DXB/TANF ( ( BETAR1+3ETAR2 ) /2 . > 

Y2A=YA-DYA $ Y2B=Y3+DYB 
TEST FOR INTERSECTION OF RAYS 
CHECK TO SEE IF CRITICAL ANGLE LIMITED. 

I COUNT =ICOUNT+l $ I F ( I COUNT . LT . 2 00 ) 244,530 

244 IF( 180.-BRGA) 245*520,255 

245 I F ( Y2A-Y2B ) 246,258,241 

246 IF( (Y2B-Y2A) .LT.YDIFF) 258,257 

255 IF ( Y2B-Y2A ) 256,258,241 

256 IFC (Y2A-Y2B) .LT.YDIFF) 258,257 v 

257 X2 = X2-XXI NC $ XX I NC= . 1*XX I NC S GO TO 241 

258 XACT = X2 $ YACT =Y2A ■ $ GO TO 400 

ROUTINE TO GET XB EQUAL TO XA WHEN ABSOLUTE VALUE OF XI IS 
LESS THAN ABSOLUTE VALUE OF XAPP 

280 IFCPHI-180. )281, 281,291 

CALCULATION OF STARTING POINTS OF RAYS WHEN PHI IS LT.' 180 

281 DX=XA-XB 

BETARl=ACOSF(VA*COSF(BETAR ) /VB) 

DY=DX/TANF ( ( BETAR+BETAR 1 ) /2. ) 

285 YB=YB+DY $ XB=XA $ V1=V0+VGRAD*XA 

ALPHA1=ALPHA $ BETA1=BETAR1/RPD 

AT THIS POINT BOTH RAYS HAVE THE SAME X COORDINATE. 

286 GO TO 300 

CALCULATION OF STARTING POINTS OF RAYS WHEN PHI IS GT. 180 
291 DX=XA-XB 

NOT£ THAT SIGN OF DX IS ALWAYS +. 

ALPHAR 1=2. *P I -ACOSF ( VB*COSF( ALPHAR )/VA) 

DY=DX/TANF ( ( ALPHAR+ALPHAR 1 ) /2. ) 

295 YA=YA-DY $ XA=XB $ V1=V0+VGRAD*X B 

ALPHA 1 = ALPHAR1/RPD $ BETA1 = BET A 

AT THIS POINT BOTH RAYS HAVE THE SAME X COORDINATE. 

296 GO TO 300 

299 V1=V0 $ ALPHAR 1 = ALPHAR $ BETAR1=BETAR 

VA=VO+VGRAD*X 1 $ VB=V0+VGRAD*XX1 

ALPHA1=ALPHA $ BET Al = BET A 

AT THIS POINT THE RAYS ARE DISTANCE STEPPED IN THE + OR -X 
DIRECTION (DEPENDING UPON PHI) UNTIL THEY INTERSECT. 
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300 DX=0.0 $ Y2A=YA $ Y2B=YB 

301 RATI OA = COSF (ALPHAR1J/V1 $ RAT I OB=COSF ( BET AR 1 ) /V 1 

XX I NC=X INC 

305 IF(PHI.LT.'180.I306»336 



CALCULATION OF XACT * YACT WHEN PHI IS LESS THAN 180 DEG. 

306 DX=DX+XXI NC 

307 V2=V1+VGRAD*DX 
IF(ABSFCV2*RATIOA) .GT.l. ) 327,308 

308 IF(ABSF(V 2* RATIOS) .GT.l.) 327,309 

309 ALPHAR2=ACOSF(V2*RATIOA) $ BFTAR2 =ACOSF ( V2*RAT I OB ) 
DYA=DX/TANF( ( ALPHAR 1+ALPHAR2 ) /2. > 

DYB=DX/TANF ( ( BETAR 1+BETAR2 ) /2 . ) 

313 Y2A = YA+DY A $ Y2B=YB+DY8 

IF(BRGA) 325*314,314 

TEST FOR INTERSECTION OF RAYS 

314 IF( 180.-BRGA) 325*520*315 

315 I F ( Y2 A-Y 2 B ) 316*328,306 

316 IF( (Y2B-Y2A) .LT.YDIFF) 328,327 

325 IFIY2B-Y2A) 326,328 ,306 

326 IF( (Y2A-Y2B) . LT . YD I FF ) 328 , 327 

327 DX = DX-XX I NC $ XX I NC = . 1*XX I NC $ GO TO 306 

328 XACT =XA+DX $ YACT =Y2A $ GO TO 400 



CALCULATION OF XACT, YACT WHEN PHI IS GREATER THAN 180 DEG. 

336 DX = DX+XX I NC 

NOTE THAT SIGN OF DX IS ALWAYS +. / 

337 V2=V1-VGRAD*DX 
IF(ABSF(V2*RATIOA) .GT.l. ) 357,338 

338 IFCABSF(V2*RATIOB) .GT.l. > 357,339 

339 ALPHAR 2=2. *P I -ACOSF( V2*RATIOA) 
BETAR2=2.*PI-ACOSF(V2*RATIOB) 

DYA = DX/T ANF ( (ALPHAR1+ALPHAR2 ) / 2 . ) 

DYB=DX/TANF ( ( BETAR 1+BETAR2 ) /2 . ) 

NOTE THAT SINCE SIGN OF DX WAS KEPT + THROUGHOUT, SIGN OF 
DYA AND DYB WILL BE THE OPPOSITE OF WHAT THEY SHOULD BE. 

343 Y2A=YA-DYA $ Y2B=YB-DYB ✓ 

TEST FOR INTERSECTION OF RAYS 

344 I F ( 180.-BRGA ) 345,520,355 

345 I F ( Y2A-Y2B ) 346,358,336 

346 IF( (Y2B-Y2A) .LT.YDIFF) 358,357 

355 I F ( Y2B-Y2A ) 356,358,336 

356 I F ( (Y2A-Y2B) .LT.YDIFF) 358,357 

357 DX = DX-XXINC $ XX I NC = . 1*XX I NC $ GO TO 336 

358 XACT =XB-DX $ YACT =Y2A S GO TO 400 

ROUTINES FOR CALCULATION OF ERRORS 



400 XERR=XACT-XAPP $ YERR=YACT-YAPP 
405 ERROR® SORTF (XERR*XERR + YERR*YERR) 
410 RACT = SQRTF (XACT*XACT + YACT*YACT ) 
415 RNGERR = RAPP - RACT 
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Ac 

c 
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c 

c 

c 

c 



SUBTRACT RNGERR FROM RAPP TO GET RACT. 

416 RELERR=ERROR/RAPP 

Z IS A DUMMY VARIABLE TO ENABLE USE OF SUBROUTINE ANGLES. 

420 Z=0.0 

421 K=3 $ CALL ANGLES < Z , Z * XACT , YACT , GAMMA , K ) 

422 DELTA = PHI - GAMMA 

SUBTRACT ERRORS FROM PHI TO GET GAMMA. 

425 K=4 $ CALL ANGLES ( Z , Z * XERR , YERR ♦ TAU , K ) 

427 BRGT = TAU - OMEGA 

428 IF(BRGT)429»500*500 

429 BRGT=BRGT+360 • 

CALCULATIONS COMPLETE* PRINT RESULTS. 

500 PRINT 510, BRGP, RACT, V2,BRGA,BRGB, ERROR, BRGT, DELTA, 
1RNGERR *RELERR 

510 FORMAT < 2X , F3 . 0 , 2X , F7 . 0 , 4X , F6 . 1 , 3X , 2 < F5 . 1 , 3X ) , F7 . 0 , 
13X,F5.1,3X,F4.1,2X,F8.0,3X,F6.4) 

515 GO TO 600 

520 PH I =PH I + 1 0 . $ PRINT 521, PHI 

521 FORMAT! F5 . 0 , 24X , 28H I MPOSS I BLE FIX(BRGP=180 DEG)) 

GO TO 600 

530 PRINT 531, PHI 

531 FORMAT < F5 . 0 , 24X , 22HCR I T I CAL' ANGLE LIMITED) 

600 CONTINUE 

INCREMENT RAPP AND GET NEW RESULTS. 

610 RAPP=RAPP+RAPPINC 

TEST TO SEE IF RAPP HAS EXCEEDED RAPPMAX. IF YES, READ 
NEW DATA, IF NOT, GET RESULTS FOR NEW RAPP. 

620 IF(RAPP-RAPPMAX) 170,170,120 
900 CONTINUE 
END 

SUBROUTINE TO CALCULATE ANGLES FROM 0 TO 360 DEGREES 
USING TANGENTS. 

SUBROUTINE ANGLES ( A , B , C , D , ANGLE , K ) 

700 RPD = .01745329252 

701 IF(D-B) 710,705,710 

705 I F ( C-A 1706,738,707 

706 ANGLE= 270 • $ GO TO 740 

707 ANGLE=90. $ GO TO 7.40 

710 ANGLER= ATANF ( ( C-A ) / ( D-B ) ) $ ANGLE=ANGLER/RPD 

715 I F ( (C-A)/(D-B) 1716,718,717 

716 IF (C-A 1720,730, 722 

717 IF(C-A) 722,736,740 

718 IF(D-B) 730,738,736 

720 ANGLE=ANGLE+360. S GO TO 740 

722 ANGLE=ANGLE+180. $ GO TO 740 

730 ANGLE= 180. $ GO TO 740 

736 ANGLE=0. $ GO TO 740 
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738 GO TO( 750*760,740,770) ,K 
740 RETURN 

750 PRINT 751 $ GO TO 740 

751 FORMAT! 35X i 18HTARGET AT SENSOR A) 

760 PRINT 761 $ GO TO 740 

761 FORMAT ( 35X, 18HTARGET AT SENSOR 8) 

770 PRINT 771 $ GO TO 740 

771 FORMAT ( 35X » 19HTARGET AT XAPP.YAPP) 
END 

END 

FINIS 

-EXECUTE. 

1500 .0001 .1 

100000 100000 500000 100000 

200000 100000 1000000 100000 

400000 100000 2000000 200000 

-999. -999. -999. -999. 



10000 

10000 

10000 

-999. 



APPENDIX II 
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Effect of Sensor Separation (SS = 10 km.) 
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Effect of Sensor Separation (SS = 100 km.) 
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Effect of Sensor Separation (SS =400 km.) 
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Effect of Target Range (RAPP = 1 km 
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Effect of Target Range (RAPP = 3 km.) 
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Effect of Target Range (RAPP = 100 km.) 
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Effect of Target Range (RAPP = 1,000 km.) 
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Effect of Target Range (RAPP = 2,000 km.) 
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Effects at Reduced Range Equal to .5 (OMEGA 
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Effects at Reduced Range Equal to .5 (OMEGA = 90°) 



APPENDIX III 



o 

m 



cc ro o r- r- m o o cnj ro ■— m m r- m o m ro 
a: *o oo O *— c\i ro ro ro cnj •— o co o cnj 

LUC or r r-f-r-f r-*— r-f— r-OC'O?' 
JOOOOOCOOOOOOOOOOO 
UJ 
O' 



romou>— in in •— ro cnj ro I s - r — oro 

cnj sr o co O *— cnj m ro ro cnj •— o co o CM 

C O ppi — r— i — »— » — •“ » — r~ • — 

ooooooooooooooooo 



II 

o 

z 

CC 

o 

LU 

cc 



II 

LL 

LL 



c^i—cNj^tr v -r-rO'or-cor-'Oroor s -^tcNjr- 

cc I I I I I J I I 

UJ I I I I I I I I I 

O 
2 : 

CC 



f— CNJ r- O ro O CO O ro r- r^- J- CNJ r— 

Jill r- f -r-r- ff -.f-f- r - r - I J | | 

I I I I I I I I I 



a 

> < 

h- r— ro^in'Or w h-czDcocor^-h — om.it ror— ro.^m'Or'-r s -cDcDCDr s -r s — om^tror- 



uj l I l I l l l l I I I l l I I l l 

o 



o 

o 

o 

o 

o 

CNJ 



II 

a. 

c_ 

< 

cc 



o 

o 

o 

o 



II 



hOOOOOOOCOOOOQOOOO CJOCOCOOOOOOCOOCCS 

— 

QfOOOOOOOOCOOOOOGOOOOOOOOOOOOCOOOOOOO 
CO o CNJ ro in O CO C7^ o CNJ l>0 ^3“ tn O UJ O r— CNJ ro LO O r- CD 0> o r— CNJ ro ^3" in nO 
r— '» t— 1 — r— r— r— r— » CNJ CNJ CNJ CNJ CNJ CNJ CNJ Qr— c— r— » r— r-i ■ r— » c— CNJ CNJ CNJ <NJ Oy{ CNJCNJ 

o 

CD 

CC O' <3 O' CNJ r- in CD CNJ h- O LO St o I s - CO r- ^ r- •— 00 1^-' >o LO O I s - CN J CO in r- CM o -=r o 

OOCNJj3TO'OCNJr— fOOCvJO^a^CMO^'— OLD || LOOr— 0>CNjCNO^CNJ>OrO»— CNJ>orOJ3*CNJ>o 

cc -d- o ro I s - o roin o o o j* cni o o ro O' a. ^j- O' m o o cnj o *o o lo ro o r- ro O' j# 

QC f-r- CNJ CNJ CNJ CNJ CNJ CNJ CNJ CNJ CNJ r-r- O ^r- CNI CNJ CNI CNJ CNJ CNI CNJ CNJ CNI 

UJ 7 CC 

CO 



a 

< 

oc 

o 

> 



o 

o 

o 

o 

o 

-=t 

II 

00 

co o 

• 

o 

o 

in 



II 



X 

LL 

CQOOOOOOOOOOOOOOOOO ooooooooooooooooo 

O • 

armomomomomomomomom-jinomomomomomomomam 
cd f— CNJCNjrorojj , ^rmmo>or w r s -cocDCDr-cDcDO^o>oO'— ^cNicNjmroj^jrmm 

•— *Cvjcnjcm CNjCNjrorororororororo rororo ro 
co 
co 
o 

< 000000000000000000-00000000000000000 

. • • 

czmomomomomomomomom*-Hmomomomomomomomom 

cQO'OO*— f-CNJCNjroro^-a-mmNOOh-r^- coocnoo*— c\JCNjroro^t-=rmm>o<3 

r— r— r— f— r— r— r-r-f— f— r— c— CNJ CNJ CNJ CM CM CNJ CNJ CNJ CNJ CNJ CNJ CNJ CNJ CNJ 



f^-CDrOCNJI^-CDNOCNJr-rOOCDO'O'OCNjrO roCNJ'OOOCOO'mi^CNJ'OCDr-CNJrOCDr- 



cjocor^-mcMO'OroG^ocMCT'j — -S'CNii — o o*— cM-zt-r-ocM'Oorooa'CNjmr-cDo 

> r -f-»-»-»- OOO OOOCDCDC0 C DCDCD CDCDCDCDCQCOO'O'O'OOC'— 



o 

o > 

O' 



II 

< 

o 

UJ 



i—r-cNij±rwr—fomh-cDN-mmoN..a-(\]*— 

0C000r-r-*-r.r-rT-,-r-0000 

<00000000000000000 

^ooooooooooooooooo 

ooooooooooooooooo 

CsJ CNJ CNJ CNI CNJ CNJ CNJ CNJ CNI CNJ CNJ CNJ CNJ CNJ CNJ CNJ CNJ 



r-oj^tf^-*nm>ON~oor^'Oror-r-^<Njr- 

OOOOr-r— r-^-r-r— r-r-r-QOOO 

ooooooooooooooooo 
ooooooooooooooooo 
oocooooaoo o o c oooo 

CNJ CNJ CNJ CNJ CNJCNJ CNJ CNJ CNJCNJ CNJ CNI CNJ CNJ CNJ CNJ CNJ 



oooooooooooooooooooooooooocooocooooo 
cc c- cnj ro m o co o o cnj po ^t m o n- o o *— cnj ro # m >o i^- cd o o •— cnj ro m 
CD r— f— •—*— f--»— f-* CNJ r-* CNJ CNJ CNJ CNI CNJ CNJ CNJ CNJCNJ cnj ro no rororo ro 



83 
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